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Abstract 


The  method  of  separation  of  variables  for  partial  differential 
equations  is  formulated  in  terms  of  direct  products  of  Hilbert  spaces. 
By  this  formulation  we  can  determine  the  possible  boundary  conditions 
for  which  the  method  is  applicable.  Finally,  a  contour  integral  repre- 
sentation 3  s  given  for  the  solution  of  some  partial  differential  equa- 
tions. 
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1.  Introduction 

The  method  of  separation  of  variables  is  one  of  the  few  methods  by 
which  partial  differential  equations  can  be  solved  explicitly.  The  method, 
however,  is  obviously  limited  to  those  cases  where  the  partial  differential 
equation  expressed  in  a  suitable  coordinate  system,  has  solutions  which  are 
products  of  functions  of  the  separate  coordinates;  consequently,  much  research 
has  been  devoted  to  the  determination  of  the  coorxiinate  systems  in  which  parti- 
cular equations,  such  as  the  wave  equation,  have  separable  solutions.  Along 
with  this,  there  has  been  a  study  of  the  special  functions  defined  by  the  sep- 
arate differential  equations.  However,  even  when  the  separable  coordinate 
systems  and  the  special  fimctions  needed  have  been  determined,  there  is  still 
one  question  to  be  discussed,  namely  the  question  of  the  type  of  boundary  con- 
ditions for  which  the  method  of  separation  of  variables  can  be  used.  To  answer 
this  question,  it  was  found  convenient  to  place  it  and  the  vrtiole  method  of  sep- 
aration of  variables  in  an  abstract  framework. 

Using  the  concept  of  the  direct  product  of  Hilbert  spaces  ,  we  show  in 
Theorem  I  that  the  idea  of  separation  of  variables  leads  to  a  method  for  solving 
operator  equations  in  which  the  operator  is  a  sum  of  direct  products  of  operators 
in  the  factor  spaces.  The  formulation  of  Theorem  I  has  many  advantages.  First, 
it  enables  us  to  discuss  how  boundary  conditions  affect  the  use  of  separation  of 
variables.  For  example,  we  show  that  certain  boundary- value  problems  for  hyper- 
bolic equations  have  a  well-defined  solution.  Also,  some  problems  vhich  do  not 
seem  to  be  amenable  to  the  method  of  separation  of  variables  can  be  solved  with 
the  aid  of  Theorem  I.  Finally,  this  theorem  suggests  a  formal  method  of  solving 
partial  differential  equations  with  a  minimum  of  calculation. 


The  concept  has  been  used  by  Cordes'-  J  in  studying  the  question  of  the 
completeness  of  the  products  of  the  separated  solutions.  However,  the 
ideas  and  methods  of  his  paper  differ  completely  from  those  of  the  present 
paper. 
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In  Theorem  II  we  consider  a  different  aspect  of  the  method  of  separa- 
tion of  variables.  In  many  problems  which  are  solved  by  this  method  there 
exist  alternative  representations  for  the  solution.  For  example,  in  solving 
Poisson's  equation  for  a  square  with  zero  boundary  conditions,  the  solution 
may  be  expressed  either  as  a  sum  of  sines  of  multiples  of  x  or  as  a  sum  of 
sines  of  multiples  of  y.  Each  representation  may  be  obtained  from  the  original 
problem  by  separating  the  variables  in  a  different  order.  In  Theorem  II  we 
obtain  a  contour-integral  representation  of  the  soiutuion  which  by  suitable 
defomnation  of  the  contour  gives  either  one  representation  or  the  other. 
Theorem  II  also  shows  that  the  Green's  function  for  a  partial  differential 
equation  may  be  written  as  the  convolution  of  the  Green's  function  for  the 
separated  ordinary  differential  equations.  Finally,  a  study  of  the  form  of 
the  solution  given  in  Theorera  II  sheds  light  on  the  role  played  by  the  Sommer- 
feld  radiation  condition  in  ensuring  the  uniqueness  of  the  solution  of  boundary- 
values  problems  for  infinite  domains. 

2.  The  abstract  formulation 

Let  ytj^^  be  a  separable  Hi Ibert  space  containing  elements  a, a-^^a^,.. ., 
u,u  ,Up,...,  f  and  let  A^,  A„,...,A^  be  self-adjoint  operators  on  /V',.   Let 
/^-  be  another  separable  Hilbert  space  containing  elements  bjb.jb^,  ,,.,v,v-,V2...  g 
and  let  B-,B_,...,B,  be  self-adjoint  operators  on  /^-p*  We  form  the  product  space 
of  /r  -  and  ;4^  in  the  manner  described  by  Murray  and  von  Neumann'--'.  Denote 
this  product  space  by  /4^,  ®  ?4^p  and  the  elements  in  it  by  w,  h,  a  ®  v,  u@b, ..., 
where,  as  stated  before,  a  and  u  belong  to  ^-t»   ^^^  "^  ^'^'^   t»  belong  to  i^^. 

The  scalar  products  in  the  spaces  /^-,   /^  and  Mi^s^p  ^^•'■^  ^®  denoted 
by  the  symbols 

<u,f>l,  <^v,g>2,   <w,h>^2 
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respectively.  The  scalar  product  in  rf ,  ©  /i^  is  related  to  the  scalar  pro- 


ducts in  A^,  and  f^^   as  follows:  If 


ra 

1 


n 


h=  z.  ^i  0gi» 


Y      J    J 


then 


m   n 

(1)  <w,h>^2  =  E  E  <v^j>i  •  <vgj>2- 

Note  that  the  scalar  product  with  the  subscript  1  is  defined  only  for 
two  elements  in  /4^  ;  however,  if  w  is  the  element  of  f4^   ®  rr    defined  above 
and  if  u  is  any  element  of  tf^  ,  we  put 

m 

(2)  <v,w>2  «  i:    <u,u.  >  V  . 

i=l 

Similarly,  if  v  is  any  element  of  /r^,  we  put 

.    .      m 

(3)  <v,w>  -  YL     <^'^i>2  \- 

j=l 

It  is  known •- J  that  if  iL.,Up,  ...form  a  complete  orthonormal  set  in 
H-  ,,  then  any  element  h  in  pf.®  /^  may  be  represented  uniquely  as  follows: 


(U) 

where 

h    -  T.  \®\. 

(5) 

\   "     <Vh>i- 

-  u  - 


Also,  the  Parsevai  equality  holds,  i.e., 

00 


(6)  <h,h>,2=  E  <\.\>2' 

Let  A  be  a  self -adjoint  operator  acting  on  a  domain  ^,    contained  in 

rT     and  let  B  be  a  self-adjoint  operator  acting  on  a  domain  Jp^   contained  in 

pfi.     We  define  the  operator  A  @  B  as  follows:  If  u^,...,u  are  elements  in 

J<,  and  if  Vt,...,v  are  elements  in  ^^,  then 
1        I'm  d 

m  m 

A0B    J2    Ui0Vi"    Ti    Au®Bv, 

We  call  A  ®  B  the  direct  product  of  the  operators  A  and  B, 

To  illustrate  these  concepts,  we  may  take  rr,  as  the  space  of  complex- 

.1 

valued  functions  f(x),  0  <  x  <  1,  such  that  J  |f(x)|   dx  <  oo,  and  take  ^- 

^»  r       2 

as  the  space  of  con^lex-vaiued  functions  g(y),   0  <  y  <  1,    such  that        |g(y)|    dy  <   oo; 

then      /r,    ®  pf-   will  be  the  space  of  complex-valued  functions  h(x,y),   0  <  x,   y  <  1, 

such  that 

[     I      |h(x,y)|^  dx  dy    <     co  . 
Jo  Jo 

The  scalar  products  will  be  defined  by  the  following  f orraxilas: 


and 


<'fp  f2>i     =       I     i7i:3n'f2(x)  dx, 
Jo 

(,&^f    g2)>2     "       I      i^^  S2(x)  dx. 


o 

1 


<[h^,    1^2^  j^2  "        f     \Cx,y)  h2(x,y)dx  dy. 
•'o 


Suppose  that 


d2 

A  u  » ^  u 

dy 
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for  ail  fiinctions  u(x)  in    Ff ^    such  that  u(x)  has  a  square-integral  second 
derivative  and  such  that  u(0)  =  u(l)  =  0,     Suppose  also  that 

d2 
B  V  •> 5"    V 

for  all  fiinctions  v(y)  in  Fr'   such  that  v(y)  has  a  square-integrable  second 
derivative  and  such  that  v(0)  =  v(l)  =  0,  Let  I  £uid  I-  denote  the  identity 
operator  in  hk  and  hk   respectively^  then,  for  all  functions  w(x,y)  in  /fi  ®  Frk 
such  that  w(x,y)  has  square-integrable  second  partial  derivatives  with  respect 
to  X  and  y,  and  such  that  w(0,y)  =  w(l,y)  =  w(x,  0)  =  w(x,l)  =  0,  we  have 


(a  ®  I^  -  I^®  b)  w(x,y)  =  -  ^  -  ^  -  -A  w. 


We  shall  now  prove  a  theorem  which  provides  an  abstract  formulation  of 
the  method  of  separation  of  variables  for  solving  partial  differential  equationSi 
Consider  a  self -ad joint  operator  L  in  >qA  ®  /^  and  suppose  that 

(7)  L=  A^  ®  B^  +  ...  +  A^  ®  Bj^ 

where  A  ,A  ,..,,A^  are  self -adjoint  operators  in  ^r^.  We  say  that  the  opera- 
tors A.^,...,A  have  a  common  spectral  representation  if  there  exists  a  non- 
decreasing  family  of  projections  E(X),  -oo  <  X  <  oo,and  a  set  of  real-valued 
functions  a  (X),,,.,a  (X)  such  that 

E(-oo)  =  0,  E(oo)  =  1, 
and  such  that  for  any  element  u  in  the  common  domain  of  A  ,.,.,A,  we  have 


(8)  AjU  «   1   a^(X)  d  E(X)u,  1<J<1^' 


-  6  - 

If  the  sf)ectrum  of  each  of  the  operators  A .  is  discrete,  the  concept  of 
a  common  spectral  representation  has  the  following  meaning:  There  exists  a  com- 
plete orthonormai  set  of  elen^ents  u,  ,u_, ...  in  /^,  and  a  set  of  real  numbers 
a.  ,  l<l<k,  l<n<  00  such  that 

(9)  ^J  ^n  =  '^jn  V  1  <  j  <  k,  1  <  n  <  oo. 

The  numbers  a   are  the  eigenvalues  and  the  elements  u  are  the  eigenelements 
of  A.,  Thus  in  this  case  a  common  spectral  representation  implies  the  existence 
of  a  common  set  of  eigenelements. 

We  now  prove  the  following: 

Theorem  I«  Suppose  that  L,  the  operator  defined  in  (?),  is  a 
closed  operator  in  rr   (x)  /tL   and  suppose  that  the  operators 
A_,,,,,A  have  a  co.Tmon  spectral  representation,  i.e.,  suppose 
(8)  holds.  Then  the  equation 

L  w  =  h, 

where  h  is  any  element  in  ^-,  ® /r^,  will  have  a  solution  if  the 
operators 

(10)        ('a^(X)B^  +  ...  +  \W\)''^ 

in  A^p  are  uniformly  bounded  for  all  X  such  that  d  E(X)  >  0, 
When  (9)  holds,  condition  (lO)  may  be  replaced  by  the  requirement 
that  the  operators 


(11)        (a^^B^.  -^^n^) 


-1 


be  uniformly  bounded  for  all  n.     If   (8)  holds,   the  SDlutuion  is  given 
by  the  formula 

(12  )  w  -     J    d  E(X)  (^(a^(X)B^  +   ...   +  \(>^)Bi^"''"  ^ 

and  if    (9)  holds,   the  solution  is 

(13)  w  =     Z    %®  (^n  ^1  "   -   *  \n^^'^  \ 

n 

where  b  is  defined  by  (5). 
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Since  the  proof  for  (i2)  involves  only  a  slight  modification  of  the 
proof  for  (l3),  we  shall  prove  (13)  only.  Consider  the  sum 


w 
mp 

n=m 


Because  the  u  are  orthonormal,  we  have 
n  ' 


y"     u    ©  (a^    B.    +   ...   +   a     B  )"^  b   . 

*—       n  '^^  \  In  1  kn  k/         n 


<%p>%>12  =  r      <t\n\  *   •••    -^  \n\)"'  ^n'      ("in^l^   -•    *  ^n^k)""^!) 

<^h     <V^n)2 

n=m 

if  M  is  the  uniform  botond  for  the  operators  in   (ll).     From   (b)  we  see  that 


00 

E/b  ,b    ^„  <   CO     , 
,        ^  n*    n'^  2 


Therefore 


^      <V^n>2 

n=m 

and  also  w       converges  to  zero.     This  implies  that 

w     =     y^     u     6^  (a.    B^    +   . . .   +  a     B,  \^  b 

m         ^       n  v>  \  In  1  kn  k/         n 

approaches  a  limit  aa  m  -*    oo.     Call  this  limit  w. 

It  must  still  be  proved  that  Lw  =  h»     We  have 

m  ^  .    ^ 

A,®  B^  %  "     r  Ai  u^  ®  B^    (a,,B^   -  ...   *  a^Bj'^  b^ 

m  J 

=     T-  "n®^n\("ln\  ^   -  "  \nV       ^n  ' 

Therefore 

L  ^  -     ^    u^®  (-i,B^  -   ...    *  -,,Bj^    (-lA  *   -   V^k)"'  ^n 

m 
=     5~    u   ®b  . 


-  8  - 

n 
Since  the  elements  u  form  a  complete  set  in  Pr~. ,  the  sim  ^  u  ®  b  con- 
verges to  h  as  m  approaches  infinity.     By  hypothesis  L  is  a  closed  operator. 
Since  w  converges  to  w  gmd  L  w  converges  to  h,  this  implies  that  L  w  ■  h. 
This  completes  the  proof  of  our  theorem. 

It  is  well  known'--'  that  if  the  operators  A..,,,., A,  commute  with  each 
other,  they  have  a  comnon  spectral  representation.  Using  this  result,  we  may 
obtain  a  generalization  of  Theorem  I, 

Let  D,,,,.,  D,  be  operators  in  rr •.   such  that  D"  has  a  bounded  inverse 

and  such  that  the  operators  D^  D-,,,,,d"  D.  are  self -adjoint,  commute,  and  have 

-1  -1 

discrete  spectra.  Then  if  we  put  Ap=  D^'  D-,.,,,  A,"  D^  D^  we  see  that  there 

exists  a  complete  set  of  orthonormal  vectors  u,  ,Up,,,.such  that 

A .  u  rt        .        ,        -, 

J         -a.     u,  2<j<k,     l<n<oo, 

''njnn*  _'J_»—  » 

and  such  that 

(lU)  D.u=a.   D^u,  2<j<k, 

^     '  jn        jnln*  — *'  — 

Using  (lU)  and  the  methods  of  Theorem  I,  we  prove  the  following: 

Corollary.  Let  D^,,.,,D,  be  the  operators  described  above  and  let 
B,,»,»,B.  be  self-adjoint  operators  in  rr  ,     Suppose  that 

M  =  D^  ®  B^  +  ..,  +  Dj^  ®  Bj^ 

is  a  closed  operator  in  M',   vy   Z^,  Suppose  that  condition  (11 ) 
is  satisfied  with  a.  defined  bv  (lU)«  Then  the  equation 

t    I 

M  w  =  h  , 

where  h  is  any  element  in  'tI  ©  ^,  has  the  solution 

where 

b  »   ^D-  u  ,  h  )>-  , 
n    \  1  n'   -^1 
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3.  Applications 

We  have  said  that  Theorem  I  is  an  abstract  formulation  of  the  method 

of  separation  of  variables.  To  see  this,  suppose  that 

^2    ^2 
L» ^ 5",  0<x,    y<l 

ex''  bj" 

on  the  domain  of  complex- valued  functions  w(x,y)  such  that  their  second 

partial  derivatives  — s-  and  — s-  are  square  integrabie  and  such  that 

dx      dy 

w(0,y)  =  w(l,y)  =  w(x,C)  =  w(x,l)  =  0  . 

Then,  as  we  have  seen  before,  we  may  put 

L  =  A  0  I^  +  I^  ®  B, 

d^  jj2 

where  A  =  -  — -^   B  ■  -  — ^  ,   and  where  I^   and  I-   are  the  identity  operators  in 

dx*^  dy'^  ^  '^ 

ph    and    /f'^y   the   spaces  of  sauarf>  integrabie  functions  of  x  and  y,   respectively, 

£iac3  A  and  L.    commute,   they  havb  a  common  spectral  representation.     In  fact, 

the  functions  u     '^  \/2     sin  nnx,   n  =  1,2,...  form  a  complete  orthonormai  set  of 

M^  t   and 

It     U       =     U    , 

In         n* 

A  2     2 

A     u    =  n     n     u  » 
n  n 

From  Theorem  I  the  solution  of  the  equation  Lw  »  h  is 

00  _  «      - 

(15)  w  -  2  7"  (n'^  n^  +  B)"-^  b  (y)  sin  nnx, 

1  "^ 


where 


b  (y)  "  I  h(x,y)  sin  nnx  dx. 


n'      I 
o 


Put 


(j?  n^  *  B)'^   b^(y)  -  c^(y). 
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Then 

(16)  2  +  n*^  n*^  c  -  b 

and  c  (0)  =  c  (1)  =  0.  The  solution  of  the  equation 

2      2 
_  d_w  _  a_w  „  j^ 

with  the  boundary  conditions 

w(0,y)  -  w(l,y)  -  w(x,0)  =  w(x,l)  •=  0 
is  found  by  the  method  of  separation  of  variables  to  be 

w  ■  2  J3  *^n^^^  ^^^  ""^' 

where  c  satisfies  (l6).  Since  this  solution  is  the  same  as  (l5 )»  we  see 
n 

that  Theorem  I  contains  the  method  of  separation  of  variables. 

Theorem  I  can  also  be  used,  however,  to  solve  certain  problems  which 
cannot  be  solved  by  separation  of  variables.  For  example,  consider  the 
following  problem: 

Find  a  function  w(x,y)  such  that 

S^w  a^w  ^K 

(17)  L  w  -  2_f  +  2  -V^r  +2_2i-h,  0<x<l,     0<y<co, 

ax^   ax'^ay^  ay^  -  -    - 

and  such  that 

w(x,0)  -  w  (x,0)  =  0  ,   w(0,y)  =  w(l,y)  -  0 


and  either 


or 


a)  w^(0,y)  =  w^(l,y)  -  0 


b)  w^(o,y)  -  W3&(i»y)  "  0- 
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Eq.  (1?)  cannot  be  solved  by  the  method  of  separation  of  variables  because, 
if  we  assume  X(.\)Y(y)  is  a  solution  of  the  homogeneous  equation,  we  get 

y(IV)     y"  »        ^(IV) 

and  clearly  in  this  equation  the  function  of  x  is  not  separable  from  the  fimc- 

tion  of  y.  However,  (l?)  can  be  solved  by  the  use  of  Theorem  I  if  the  boundary 

d^  d^ 

conditions  b)  are  assumed.  To  see  this,  p\rt  A  ■  — p  and  A-  *  -  — ^  with  the 

■^   dx^     '^     dx"^ 

boundary  conditions 

u(0)  =  u(l)  '   u"(0)  »  u"(l)  «  0 

and  put  B-  »  -  — j  ,  B,  ■  — r-  ,  0  <  y  <  oo  with  the  boimdary  conditions 
dy'^    ^       dy^ 

v(0)  "  V  (0)  =  0.  Then 

L  -  A^  ©  Ig  +  2  Ag  ®  B2  +  1^0  B^. 

The  operators  A^,  A-  and  I^  have  a  common  spectral  representation  defined  by 
the  complete  set  of  eigenf  unctions  u  »  \/2"  sin  nnx,  n  »  1, 2,  • . .  •  Since 

2  2 
A_u  »  n  n  u 
2  n        n 


and 


A-u  «  n  n  u  , 
in        n* 


we  need  only  show,  in  order  to  be  able  to  apply  Theorem  I,  that  the  operators 

(B3  *  2  n\\  .   n^n^)-l 

are  uniformly  bounded.  The  proof  of  the  existence  of  t^iis  boxuid  follows  easily 
from  the  fact  that  the  operators  Bp  and  B-  are  positive-definite.  Using  (l3) 
we  may  obtain  an  explicit  represen-t  ation  of  the  solution  for  (17).  The  result 
is  not  necesssry  for  our  purposes;  the  interest  of  this  problem  is  that  it 
cannot  be  solved  by  separation  of  variables  but  can  be  solved  by  Theorem  I» 


-  12  - 


Note,  however,  that  if  we  assume  boundary  conditions  a).  Theorem  I 

d^ 
cannot  be  used.  We  may  define  B^  and  B-  as  before,  but  now  A.^  =  — j-  and 

A_  ■  — ^  f   with  the  boundary  conditions 
'^       dx 

u(0)  =  u(l)  =  u  (0)  «  u(l)  =  0. 

These  operators  A.^  and  A^   do  not  have  a  common  spectral  representation  and 
therefore  our  methods  do  not  apply.  In  fact,  an  explicit  representation  of 
the  solution  of  (i?)  with  the  boixndary  conditions  a)   is  not  known. 

Another  important  consequence  of  Theorem  I  is  that  it  provides  a  proof 
of  the  existence  of  a  solution  to  the  equation 

L  w  »  h. 

However,  this  existence  is  proved  only  if  a  condition  such  as  (10)  or  (11)  is 

satisfied.  It  is  interesting  that  tnis  condition  may  be  satisfied  even  in 

cases  where  the  general  theory  of  partial  differential  equations  seems  to  state 

that  no  solution  is  possible.  For  example,  according  to  the  general  theory,  a 

boundary-value  problem  for  a  hyperbolic  equation  is  an  unreasonable  problem. 

However,  consider  the  following: 

Find  a  function  w(x,y)  satisfying  the  hyperbolic  equation 

2     2 
(18)  Lw=  ^-  ^^=  h(x,y) 

Sx"   dy^ 

and  satisfying  the  boundary  conditions 

w(0,y)  =  w(i,y)  =  w(x,C)  =  w  (x,l)  =  0. 

If  we  put 

L  =  A  0  I^  +  I^  0  B, 

^2  ^2 

where  A  = r-  ,  with  the  boundary  conditions  u(0)  =  u(l)  =  0,  and  B  «  +  — sy- 

dx*^  ,  dy"^ 

Tvith  the  boundary  conditions  v(0)  =  v  (l)  =  0,  then  we  find  that  the  eigenvalues 

2  2  2  2  ""l 

of  A  are  n  it  and  the  operators  (B  +  n  n  )   are  uniformly  bounded.  Consequently, 
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Theorem  I  may  be  applied  and  Eq.  (l8)  has  a  solution.  However,  if  we  change 
the  boundary  condition  w  (x, 1)  =  0  to  w  (x,l  +  e)  =  0,  then  there  exist  arbit- 
rarily  small  values  of  e  for  which  Eq.  (l?)  will  not  have  a  solution.  This 
type  of  behavior,  in  which  an  infinitesmai  variation  in  a  parameter,  such  as 
e,  may  cause  a  problem  to  change  from  one  having  a  solution  to  one  not  having 
a  solution  and  then  back  again,  is  physically  unreasonable  because  no  parameter 
can  ever  be  measured  exactly.  This  is  the  basis  for  the  theoretical  conclusion 
that  the  boundary- value  problem  for  a  hyperbolic  problem  is  unreasonable. 

Theorem  I  has  an  important  practical  consequence.  The  form  of  the 
solution  for  w  shows  that  the  operators  a  B  +  ...  +  a  B,  must  be  inverted. 
This  means  that  we  may  consider  the  operators  A  ,  ,.,,  A,  as  constants,  invert 
the  operator  A.^B^  +  ...  ■••  A.  B,  (here  A^,...,  A,  represent  constants),  and  then 
interpret  the  result  by  using  the  spectral  representation  for  A.^,...,A,  , 

An  illustration  will  clarify  this  method.  Consider  the  problem  of  find- 
ing a  function  w(x,  y)  such  that 


(19) 


and  such  that 


2     2 
d  w   d  w   ,  /    X 
9 9  "  h(x,y) 


0  <  X  <  1,  -00  <  y  <  00 


w(0,y)  »  w(l,y)  =  0, 

d2 
Ptit  A  -  -  TT  ,  with  the  boundary  conditions  u(0)  -  U(l)  =  0.  Then 

(19)  may  be  written  as 


2 
d_w 


-  A  w  »  -  h(x,y)  . 


The  solution  of  this  equation  (assuming  w  is  not  exponentially  large  at  y  =  +  oo)  is 


'-00 


-  lU  - 


In  order  to  interpret  this  result  we  may  use  the  operational  calculus  of 
Dunf ord  L-^-'  or,  what  is  equivalent,  the  spectral  representation  of  A,  The 


2_2 
runcxions   oi    a  are   u     =  \jc      sin  nnx  cina   bne   corresponaing  are  r 

Since 


eigenfunctions  of  A  are  u  =  s/^    sin  nnx  and  the  corresponding  are  n  n  . 


f (A)u^  =  f (n^n^)  u^. 


and  since 

00 


h(x,  >7)  =  2  ^  sin  nnx     f     sin  nnx  ?  h(^,  /7  )d4, 

we  see  that 

1  00 

w  =  r    2in-2H     r    sin  nn4  d4     J       e"""!^"^  lh(C,  //  )d  /J  . 

Thus  the  desired  solution  to   (19)  has  been  obtained  with  a  minimum  of  calcu- 
lation. 


U,     The  Euler-Poisson-Darboux  equation 

The  methods  of  the  preceding  section  can  be  applied  to  the  discussion 
of  the  Euler-Poisson-Darboux  equation'--',  namely, 

(20)  A  w  =  w^^  +  k  t"-"-  w^, 

with  the  initial  conditions  w  =  f  and  w.  =  0  for  t  =  0.  Here  the  Lapiacian 
is  an  operator  in  m  space  variables.  For  convenience,  put  A«  -p  ,  and  con- 
si  ier  :vq.  (20)  with  p  a  constant.  It  becomes 

(21)  w^^  +  k  t"-*-  w^  +  p^w  »  0. 
The  solution  of  tliis  equation  is 

w  =  a^t'^  J^(pt)  +  a2t'^J_^(pt), 

where 

a  «  (i  -  k)/2 
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and  where  a.  and  a  are  independent  of  t,  and  J  and  J_  are  Bessei  functions. 

To  satisfy  the  imtiai  conditions  we  must  find  a  solution  of  (21)  such  that 

w(0)  =  f  and  w. (O)  »  0.  From  the  power  series  expansion  for  the  Bessei  func- 
t 


tions  V.  have 


^>^'-7d^    (fl[-0(t^)]. 


rc'^+i) 

Consequently,  for  ^  >  j  the  function  t  J^(pt)  and  its  first  derivative  both 

1  tt  /   \ 

vanish  for  t  =  0,  For  a  <  —  ,  the  derivative  of  the  function  t  J  (pt)  becomes 

infinite  for  t  =  0.  On  the  other  hand,  the  value  of  the  function  t  J__(pt)  at 

1    /2f 
t  •  0  is     '     I —I  while  its  derivative  vanishes  at  t  ■»  0.   Hsing  these 

n(i-a)  \p/ 

results,  we  see  thit  the  solution  of  (21)  which  satisfies  the  initial  conditions  is 

(22,      ^-r^)  (fY  ^.jpt)f. 

To  interpret  this  we  must  use  the  spectral  representation  for  ZA  ■  -p  ,  The 
spectral  representation  will  be  obtained  by  using  the  Fourier  transform  in  m 
dimensions.  Let  x,  ^  denote  vectors  in  m-dimensionai  Euclidean  space  and  £ 
a  vector  of  m  dimensions  in  the  transformed  space.  Then  we  have 

f(x)  .  (2n)-"^  f...  re^2-/^-^^f(Od£  d^ 

where  the  integration  is  extended  over  the  entire  p  and  ^  space.  Also,  if 
5f(-A)  is  an  analytic  ftinction  of  ^,  we  have 

0(-A)f(x)  .  (2n)-'"  r...  r0(p2)e^e-(x-i)f(^)d£dE, 

2 

where  p  =  £  *  £•  Using  this,  we  see  that  (22)  becomes 

(23)  w  -  C^  J...J{ptfj_^{pt}   e^E'^^-^^  f  (i)di  dp. 


where 


^  .  p  (  ilil  )  2-(2„)-"  . 
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By  interchanging  the  order  of  integration,  (23)  can  be  reduced  to  the  fommias 
given  elsewhere'--'.  Simiiar  methods  have  been  proposed  in  an  abstract  by  P.O. 
Rosenbium '--*-'  , 

The  above  discussion  has  shown  that  the  Euler-Poisson-Darboux  equation 
has  a  unique  solution  for  k  >  0  but  not  for  k  <  0.  However,  for  k  <  0  we  may- 
consider  the  solution  (23)  which  is  the  analytic  continuation  of  the  solution 
for  k  >  0,  As  can  be  seen  from  the  formula  for  C,  ,  this  solution  breaks  down 
when  k  =  -1,  -3,  -5>  • • •  •  This  breakdown  could  have  been  foreseen  because 
the  exceptional  values  of  k  correspond  to  integral  values  of  a.   and  for  these 

values  of  a  the  fvinctions  J  and  J   are  not  linearly  independent. 

a     —a 

Let  us  consider  the  exceptional  values  of  k.  They  occur  when  k  =  l-2n 
or  a  *  -n,  with  n  a  positive  integer.  For  these  values  of  k  we  take  the  follow- 
ing solution  of  (21): 

(2U)  w L_  (^f  N  (pt)f, 

where  N  (pt)  is  the  Neumann  fuiiction.  This  solution  has  the  disadvantage 
that  its  2n-th  derivative  becomes  logarithmically  infinite  for  t  ■  0,  However, 
for  certain  kinds  of  functions  f (x),  the  infinity  will  not  occur;  consequently, 
w  and  all  of  its  derivatives  will  be  continuous  at  t  =  0.  To  find  these  func- 
tions, note  that 

(pt)X(pt)  '    n  C,(pt)^^  *  (pt)^''  log(pt)  2  c'(pt)2J, 

where  C.  and  C.  are  scaiars.  The  solution  w  will  not  have  any  logarithmic  term  if 

(pt)2"  f (x)  -  0, 
that  is,  if 

(25)  a"  f  (x)  .  0. 
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Thus,  we  have  obtained  Weinstein's  result"- J  that  if  f(x)  satisfies  (25),  then 
there  exists  a  solution  of  (20)  wMch  is  continuous  and  has  continuous  deriva- 
tives of  all  orders. 

5.  Alternative  representations 

Formula  (l5)  Is  that  solution  of  the  equation 

2     2 

3x    dy 

which  satisfies  zero  boundary  conditions  on  the  sides  of  the  square  x  ■  0,  x  -  1, 

y  =  0,  y  =  1.  This  solution  is  expressed  in  terms  of  the  x-operator  eigenfunc- 

tions  sin  nnx.  Since  the  problem  is  symmetric,  it  is  to  be  expected  that  there 

will  exist  a  representation  of  the  solution  in  terms  of  the  y-operator  eigenfunc- 

tions,  namely  sin  nny.  This  expectation  is  correct,  and  the  desired  representa- 

^2 
tion  can  be  foimd  by  using  Theorem  I,  since  the^  operators  -  — r-  commute  with  I^. 

'  dy*^  '^ 

It  would  be  desirable  to  obtain  both  representations  without  having  to 

return  to  the  original  problan.  This  will  be  achieved  by  Theoran  II,  in  which 

we  shall  obtain  a  contour  integral  representation  for  the  solution.  Deforming 

the  contour  one  way  will  give  a  representation  in  terms  of  x-eigenfunctions; 

deforming  the  contour  in  another  way  will  give  a  representation  in  terms  of 

y-eigenf unctions . 

Theorem  II«  Suppose  that  the  closed  operator  L  in  ^^  (E) '^o   ■"■^  given 
by  the  f  oraula 

L  =  A  (^  I^  +  I-^  0  B, 

where  A  and  B  are  self -adjoint  operators  in  /r,  and  Ar_  respectively. 
Suppose  that  the  distance  between  the  spectrum  of  -B  and  the  spectrum 
of  A  is  greater  than  zero.  Then  the  solution  of 


L  w  =  h, 

where  h  belongs  to  /t   CV  'r ^,  is 
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(27)        *  "  5^   I   (A  -  X)"^  0  (B  +  X)--^  dX  h, 

where  the  region  bounded  by  C  contains  the  complete  spectrum  of 
A  and  is  at  a  non-zero  distance  5  from  the  spectrum  of  -B. 

As  in  Theorem  I,  we  shall  restrict  the  proof  to  the  case  where  the  spec- 
trum of  A  is  discrete  and  contains  the  points  X^ ,  Xp,...  .  Let  C  be  any  closed 
contour  containing  the  points  X-,  X-, ».«,  X  and  let  the  domain  bounded  by  C 
be  at  a  distance  greater  then  5  from  the  spectrum  of  -3.  Consider  the  integral 

w  -  (2ni)"-'-  I   (A  -  X)"-"-  0(B  +  X)"-"-  dX  h. 

From  (U)  we  have 

h  =  T"  u  ®  b  , 
^—   n  ^-^  n' 

and  therefore 

w  «  (ani)""^  I   (X  -  X)"-'-u  0(B  +  X)"-H)  dX. 
p^       I   ^n       n^-'^        n 

P 

By  the  definition  of  C  ,  the  operator  B  +  X  is  bounded  by  5.  It  is  easy  to 

show  that  we  may  interchange  integration  and  summation;  then  we  obtain 

Wp  -  ^  u^  0  (B  -  X^)-'  b^. 

But  we  have  shDwn  in  the  proof  of  Theorem  I  that  w  converges  to  a  limit  w 
which  is  the  solution  of 

L  w  =  h. 
This  proves  Theorem  II, 

In  the  proof  of  Theorem  II  we  have  shown  how  (28)  leads  to  an  expansion 
in  terms  of  A-sigenf unctions.  Similarly,  if  the  operator  (A-X)"  ®  (B  +  X) 


-  19  - 


is  sufficiently  small  for  large  |X|,  we  may  deform  the  contour  C  to  enclose 

the  spectrum  of  -B  and  thus  obtain  an  expansion  in  terms  of  B-eigenf unctions. 

As  an  illustration,  consider  again  (26).  Here  A  =  -  — s-  ,  with  the  boundary 

^2  dx'^ 

conditions  u(0)  =  u(l)  =  0,  and  B  =  -  — r->  with  the  boundary  conditions 

v(0)  =  v(l)  =  0,  The  operators  (A-X)   and  (B+X)~  are  integial  operators  -.lith 
kernels  g(x,4,X)  and  g(y,   ,-X)  respectively,  where 


g(x,5,X)  .  slny^xsinVT  (j-^) 

\A  sin  yr 


X  <  4, 


sin  \/T  g  ^in  -/^   (l-x) 
y/T.  sin  V^ 


X  >  4. 


From  Theorem  II  we  have 


(28) 


""  2ii  J  "^^  j   f  g(x,S,X;g(y,  )7,-X)h(c,r|)d4  dKj, 


'C    'o  "0 

where  C  is  the  contour  indicated  in  Fig,  1, 


X- plane 


c 

c 

/""^ 

-4 

-[^y 

V^^ 

Figure  1 
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If  we  evaluate  the  residues  inside  C,  we  get  (l^).  However,  if  we  deform  C 
into  C  ,  there  will  be  poles  at  X  =  -1,  -U,  ...  and  by  evaluating  the  residues 
at  these  poles  we  would  get  a  representation  for  w  in  terms  of  sin  nny. 
Theorem  II  may  also  be  written  as  follows: 

/a  ©I2  +  1^0  bT'   -   ^   r  (A-X)-^  0  (B+X)-^  dX. 

For  differential  operators  the  inverse  operator  is  an  integral  operator. 
Let  G(n)  be  the  kernel  of  the  integral  operator  which  inverts  L  -  ^ji  and  let 
g^ (X),  g^(X)  be  the  kernels  of  the  integral  operators  which  invert  A-X  and 
B-X  respectively.  Then  Theorem  II  shows  that 

r 


G(^)  -  (2ni)"^ 


g^(X)g2(ti  -  X)  dX. 


If  we  put  h  =  5  (4-4  )5()0-  Y\   )  in  (28),  it  will  be  an  illustration  of  this 
result. 

The  following  Corollary  to  Theorem  II  may  be  proved  by  a  method  similar 
to  that  used  to  prove  the  Corollary  of  Theorem  I, 

Corollary.  Suppose  that 

L  -  A^  0  B^  +  A^  ©  B2, 

where  A^,  A  are  operators  in  H^-y   and  R,,  B_  are  operators  in   //^. 
Suppose  that  A,  and  B^  have  bounjled  inverses,  that  A,  A_  and  B-  B, 
are  self -adjoint,  and  that  the  distance  from  the  spectnun  of  A  A  to 
the  spectrum  of  B~  B^  is  positive.  Then  the  solution  of 

L  w  =  h 
is 


w  =  (2ni)"^  J   (A^-  XA^)"^  0  (B^+  XB^)""^  dX  h. 
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6.  Interlacing  spectra 

In  this  section  we  shall  discuss  the  possible  relationships  between 
the  spectra  of  A  and  -B  which  permit  the  use  of  Theoran  II.  The  simplest  re- 
lationship is  that  illustrated  by  the  problem  (26 j,  in  which  the  spectrian  of  A 
is  a  set  of  discrete  points  on  the  positive  real  axis  and  the  spectrum  of  -D 
is  a  set  of  discrete  points  on  the  negative  real  axis.  This  relationship  will 
occur  if  A  and  B  are  both  positive-definite  operators. 

If  A  is  positive-definite  and  B  is  negative-definite,  the  spectnai  of 

both  A  and  -B  will  be  on  the  positive  real  axis.  This  occurs  in  (l8),  ^ere 

2 

the  spectrum  of  A  consists  of  the  points  X  »  (nx)  and  the  spectrum  of  -B 

12  2 
consists  of  the  points  \i    *  (n  +  j)  n  ,  Since  the  distance  from  the  set  of 

X  to  the  set  of  ^    is  greater  than  zero.  Theorem  II  may  be  applied  and  will 

give  a  result  equivalent  to  that  obtained  previously.  An  appropriate  contour 

C  is  illustrated  in  Fig.  2. 


i 

^^ 

'"( 

\Mi  y 

'        Xg      ^ 

K^^y 

X3     V 

FJCTre  2 
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Let  us  see  what  happens  to  the  spectrum  of  B  if  we  change  the  length 
of  the  veticai  side  to  1  +  e  so  that  the  boundary  conditions  are  v(0)  »  v  (1+e)  =  0, 

Me   find  that  the  spectrin  of  -B  consists  of  the  points  \i.    =  '^  i  i..^  i  .  Note 
that  for  arbitrary  small  values  of  e  there  exist  values  of  n  and  m  such  that 
|X  -ti  I  can  be  made  arbitrarily  small.  This  shows  that  the  spectra  are  so 
closely  interlaced  that  we  cannot  find  a  contour  C  and  thus  we  cannot  apply 
Theorem  II.  In  fact,  for  these  values  of  e  the  problem  (16)  does  not  have  a 
solution. 

If  A  is  positive-definite  and  B  is  onlv  bounded  below,  the  spectrum  of 
-3  may  partly  overlap  the  spectrum  of  A.  Suppose  that  the  spectrum  of  A  is 
discrete,  with  its  only  limit  point  at  infinity.  If  the  spectrum  of  -B  is  also 
discrete,  with  its  only  limit  point  at  -oo,  then  there  are  two  possibilities: 
either  the  spectra  have  a  positive  distance  between  them,  or  there  is  a  value 
of  \  wliich  is  both  in  the  spectrimi  of  A  and  in  the  spectrum  of  -B.  If  first 
possibility  occurs.  Theorem  II  may  be  applied.  If  the  second  possibility  occurs, 
the  equation 

L  w  =  0 

will  have  a  non-trivial  solution,  namely  w  =  u  ®  v  ,  where  u  ,  v^  are  the 

n    n'       n'  n 

eigenelements  of  the  operators  A  and  -B  respectively,  corresponding  to  the  common 
eigenvalue  X  .  The  question  of  the  existence  of  a  solution  to  L  w  =  h  is  now 
more  difficult  and  will  not  be  discussed. 

If  the  spectrum  of  -B  is  continuous,  it  may  overlap  a  point  in  the 
spectrtm  of  A  as  is  shown  in  Fig.  3  (the  shaded  pairt  of  the  real  axis  repre- 
sents the  continuous  spectrum  of  -B). 
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Figure  3 

Since  X^  is  in  the  spectrum  of  -B,  there  is  difficulty  in  obtaining  a 
suitable  contour  C,  If  we  consider  an  ex?jnple  in  which  this  occurs,  such  as 
(19),  we  find  that  because  of  this  overlap  the  problem  is  not  wel]  posed  unless 
we  add  a  condition  at  infinity,  namely  the  Sommerfeld  radiation  condition.  IVe 
may  say  that  the  role  of  this  condition  is  to  determine  the  behavior  of  the 
curve  C  in  the  neighborhood  of  >l.  (see  Fig.  U). 
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Figure  h 
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